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ABSTRACT 

We examine Milky Way-Magellanic Cloud systems selected from the Millennium-II 
Simulation in order to place the orbits of the Magellanic Clouds in a cosmological 
context. Our analysis shows that satellites massive enough to be LMC analogs are 
typically accreted at late times. Moreover, those that are accreted at early times and 
survive to the present have orbital properties that are discrepant with those observed 
for the LMC. The high velocity of the LMC, coupled with the dearth of unbound orbits 
seen in the simulation, argues that the mass of the MW's halo is unlikely to be less 
than 2 x 10 12 M Q . This conclusion is further supported by statistics of halos hosting 
satellites with masses, velocities, and separations comparable to those of the LMC. We 
further show that: (1) LMC and SMC-mass objects are not particularly uncommon 
in MW-mass halos; (2) the apparently high angular momentum of the LMC is not 
cosmologically unusual; and (3) it is rare for a MW halo to host a LMC-SMC binary 
system at z = 0, but high speed binary pairs accreted at late times are possible. Based 
on these results, we conclude that the LMC was accreted within the past four Gyr 
and is currently making its first pericentric passage about the MW. 

Key words: Galaxy: fundamental parameters - Galaxy: formation - galaxies: for- 
mation - galaxies: kinematics and dynamics - Magellanic Clouds 



1 INTRODUCTION 

The Milky Way and its satellite galaxies offer a unique lab- 
oratory for near-field cosmology. The proximity of Milky 
Way (MW) satellites means that their stellar content is re- 
solved and can be used as a probe of galaxy formation and 
evolution ( Grebel||2005 1. Furthermore, with the high astro- 
metric precision of instruments such as the Advanced Cam- 
era for Surveys on the Hubble Space Telescope (HST), it 
is now also possible to measure accurate proper motions 
for some of these satellites (e.g., Kallivayalil ct al. 20 06a|b 



Piatek et al.||200~ 2008). These measurements have signif- 
icantly improved constraints on the satellites' orbital his- 
tories, and have also revealed some surprises. In particu- 
lar, [Kallivayalil et ah] ( |2006a| hereafter K06) found that the 
velocity of the Large Magellanic Cloud (LMC) is approxi- 
mately 380 kms -1 , which is much large r than what typically 



had been assumed (<300kms ; e.g., Gardiner & Noguchi 



1996[ ) in modeling the orbit of the LMC - 

This large velocity has forced a reconsideration of the 
conventional picture of the orbital history of the Magellanic 
Clouds (MCs), wherein the MCs have made multiple pas- 
sages about the MW over a Hubble time. K06 and |Besla| 



et al. (2007 hereafter, B07) analyzed possible LMC orbits 
and divided them into two categories: early accretion, in 
which the LMC has made at least one complete orbit about 
the Milky Way; and late accretion, where the LMC is cur- 
rently making its first pericentric passage. Distinguishing be- 
tween these two orbital histories has important consequences 
for our understanding of the Local Group's assembly history 
and for the formation of the Magellanic Stream ( |Besla et al.| 



2010) 



The high velocities of the Clouds also raise the question 
of whether the MCs' orbits are typical of massive satellites 
in MW-like systems and whether they can provide informa- 
tion about the mass of the MW's halo. B07 showed that 
the MCs are effectively on an unbound orbit if the MW's 
mass is on the low end of current estimates (~ 1O 12 M0, 
consistent with the results of IXue et al. 2008] based on blue 
horizontal- branch stars). On the other hand, a massive MW 
halo (~ [2 — 2.5] x 10 12 M©) - in line with estimates based on 
the timing argument and satellite kinematics including data 
from Leo I ( |Li fc White|2008| [Watkins et al.|20To| ) - implies 
that the LMC has a velocity that is substantially lower than 
the local escape speed. 

Analyses of cosmological simulations indicate that un- 
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2006| |Diemand et al.|2007||Wetzel|2010[ ): for example, |Wet 
zel ( 2010 1 found that less than 2% of merging satellites have 



formally unbound orbits at all masses and redshifts. |Sales| 
et al. (20071 have shown that a non-negligible number of 



subhalos in cosmological simulations can be scattered to 
high energy orbits as the result of three-body encounters, 
but these subhalos are dynamically required to be low-mass, 
likely rendering this mechanism irrelevant for the LMC. Al- 
most all of these results are based on modeling dark matter 
halos as point mass, Kepler potentials. This is a poor ap- 
proximation to the true structure of dark matter halos at 
small radii, however. The LMC, which is currently at < 
one-fifth the MW's virial radius, is in precisely this regime. 
As a result, it is far from clear whether an unbound LMC 
would be highly unusual in a cosmological context. 

The very existence of massive satellites such as the MCs 
can also be used to place constraints on the mass of the MW. 
If the LMC is typical for galaxies of its stellar mass, it likely 
had a dark matter mass of [1 — 2] x 10 11 Mq before accretion 
by the MW ( |Guo et al.||2010| |Boylan-Kolchin et al.|[2010 l. 
The SMC has a stellar mass that is typical for halos that are 
only a factor of 2-3 less massive than this. |Boylan-Kolchin| 
et al. (2010 hereafter BK10) used statistics of subhalos from 



a large sample of simulated MW-mass halos to conclude that 
MC-mass galaxies are rare if the MW has M < 10 12 M© but 
are much more typical if the MW has a massive dark matter 
halo (~ 2.5 x 10 12 ). These results were based purely on the 
masses of the MCs and MW, however, and did not include 
the additional constraints provided by the MCs' orbits. 

In this paper, we examine the orbits, accretion epochs, 
and masses of MW-MC systems using the Millennium-II 



Simulation (MS-II; Boylan-Kolchin et al. 20091, which has 
sufficient mass resolution to probe SMC-mass scales and a 
large enough volume to contain a large, statistical sample 
of MW-mass halos. We use realistic - but still spherically- 
symmetric - potentials that are expected for LCDM dark 
matter halos when computing orbital parameters in an at- 
tempt to accurately describe the orbits of MC analogs at 
small radii from their hosts. 

Our goal is to investigate five main questions related 
to the MCs in order to place their orbital properties and 
masses in a cosmological context: 

• How common are satellites with masses similar to those 
of the LMC and SMC within Milky Way-mass halos at red- 
shift zero? 

• How typical is the LMC in terms of its orbital properties 
(e.g., energy and angular momentum)? 

• Are LMC-type satellites at z = typically accreted 
early (having completed at least one pericentric passage 
about their hosts) or late (being on their first infall towards 
their hosts)? 

• How likely is it that the MCs were accreted as a binary 
system? 

• What can we infer about host halos from the properties 
of massive satellites such as the LMC? 

Our work is structured as follows. §|2] provides relevant 
details about the MS-II, describes our assumptions about 
the MW and MCs, and defines our subhalo samples. In § |3.1| 
we determine the most likely accretion epoch for the LMC 
based on mass considerations. In § |3.2f § [374] we fold the or- 
bital properties of the LMC into the analysis to distinguish 



6.2|. We present our 



between the early or late accretion scenarios. Section [1] ex- 
plores the expected frequency of L/SMC-type companions 
about MW-type hosts based on the MCs' expected infall 
masses. In § [5] we examine potential L/SMC binary sys- 
tems and their properties. Our results are discussed in § [6] 
in particular, we assess the likelihood that the MCs are on 
their first passage about the MW (§ 
conclusions in § [7] 



2 METHODOLOGY 

2.1 The Millennium-II Simulation 

Our analysis of objects similar to the Magellanic Clouds 
is based on the Millennium-II Simulation, a very large TV- 
body simulation that follows the evolution of over ten billion 
particles in a periodic cube of (137 Mpc) 3 from redshift 127 
to 0. The cosmological parameters used in the MS-II are 
identical to those adopted for the Millennium Simulation 
( Springel et al. 2005 1 and the Aquarius simulations ( |Springei] 
et al.|2008[ )~ 



fitot =i.o, n m =0.25, n b = 0.045, n A = 0.75, 

h = 0.73, (j 8 = 0.9, n s = l, (1) 

where h is the Hubble constant at redshift zero in units of 
100 km s" 1 Mpc" 1 , as is the rms amplitude of linear mass 
fluctuations in 8 ft" 1 Mpc spheres at z = 0, and n s is the 
spectral index of the primordial power spectrum. The MS- 
II offers a unique combination of mass resolution and large 
volume for studying dynamics of Magellanic Cloud analogs 
within Milky Way-mass systems: the MS-II particle mass 
of m p = 9.43 x 10 6 M e results in over 100,000 particles in 
Milky Way-mass halos at z = 0, and resolves LMC-mass 
subhalos with > 10, 000 particles. For further information 
about the MS-II, see |Boylan-Kolchin etH] ( |2009[ |p1 



2.2 Milky Way-mass halos 

Current estimates of the mass of t he Milky Way's halo 
range from « [1 - 3] x 10 12 M Q (e.g., |Sakamoto et al.|2003] 
Battaglia et al.|2005|[Dehnen et al.|2006||Xue et al.|2008| |Li| 
fc White|2008| fCnedin et al.|2010[[watkins et al.|2010| and 



references therein). In order to bracket this range, and to 
understand any trends with the mass of the MW, we select 
all halos with 4.3 x 10 11 < Af vir /M < 4.3 x 10 12 from the 
MS-II at redshift zero as our primary "Milky Way" sam- 
ple. This set is identical to that of Boylan-Kolchin ct al. 
(20101 and contains approximately 7600 dark matter halos, 
2658 of which have M vil <E [1 - 3] x 10 12 Mq. It is impor- 
tant to note that not all of these halos are expected to host 
MW-like galaxies in the standard ACDM model: for exmple, 
approximately 30% of halos in our full MW sample should 
host galaxies that are not late-type based on their colors 



and specific star formation rates (Weinmann et al. 2006) 



By taking subsets of our main sample that are based on, 



Merger trees, along with halo and subhalo catalogs, 
from the MS-II are publicly available at |http://www.mpa-| 
garching.mpg.de/galform/millennium-II/ . 
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e.g., environment, we can explore whether host halo proper- 
ties correlate with the likelihood of hosting satellite galaxies 
like the MCs. 

We use the radius inside of which the average density 
exceeds the critical density of the universe by a factor A v h- 
(see Eke et aL||1996 and |Bryan fc Norman|1998| for details) 
as the virial radius 7? v i r of our halos, and the mass enclosed 
within this radius, M v j r , as the virial mass. At redshift zero 
in the cosmology of the Millennium-II Simulation, A v i r = 
94.2 and the virial radius 7? v i r and virial (circular) velocity 
V v ir scale with the virial mass as: 



iivir = 257.8 
Kir = 129.2 



M vir 

10 12 M Q 
M. 



1/3 



kpc 



1/3 

vir \ n —1 

km s 



10 12 M Q 



(2) 
(3) 



We will use these relations in subsequent sections to scale 
measured properties of the Magellanic Clouds - e.g., their 
angular momenta - to the virial values for a range of virial 
masses for the MW. 

Later sections will also focus on orbital properties of 
subhalos within their hosts. These properties will be com- 
puted assuming that the host dark matter halos are well- 
fitted by Navarro, Frenk, and White (119961 119971 hereafter, 



NFW) density profiles. The structure of the NFW profile 
is fully specified by a halo's mass and the radius at which 
the circular velocity curve peaks, i? ma x. Using values of M v i r 
and i? max computed directly from the MS-II to set the NFW 
potential, we calculate the energy, angular momentum, and 
eccentricity of orbiting subhalos within MW-mass host^] 
We defer an analysis of the current distance of the LMC 
from the MW to Section |6.3[ as it is inherently extremely 
rare for a satellite to be near pericenter (as is the case for 
LMC) for the eccentric orbits typical of ACDM satellites 



(Tormen 1997 Diemand et al. 20071. Note that one limi- 



tation of our approach is the use of spherically symmetric 
profiles in our calculations, whereas halos from cosmologi- 
cal dark matter simulations are typically prolate or triaxial 
(e.g., [Warren et al.||1992| |Jing fc Suto||2002| [Allgood et al. 



2006 



Bett et al.|2007j 7 



2.3 The Magellanic Clouds 

2.3.1 Observed properties 

The Large (Small) Magellanic Cloud is the most (second- 
most) luminous satellite of the Milky Way. Following Kim 
et al. ( 1998 1, we assume that the stellar mass of the LMC is 
M L mc,* = 2.5 x 10 9 Mq. We adopt a galactocentric distance 
of -Rlmc = 50.1 kpc ( Freedman et al."||2001 1. K06's analysis 



of the LMC's proper motion shows that 



VLMC.ta 
VLMC.ra 



367 ± 18 kms" 
89±4kms" ] 



(4) 



2 We assume a fixed concentration of 10 when computing the 
potential of the actual MW. This is not a strong assumption, 
as changing the concentration between, e.g., 8 & 12 does not 
significantly affect the inferred pericenter, angular momentum, 
or orbital energy of the subhalos. 



The specific angular momentum of the LMC, normalized by 
the specific angular momentum of a circular orbit at the 
MWs virial radius, is therefore 



JLMC 



-Rlmc 14mc,i 

Rvir ^vir 



= 0.55 



Mvir 



10 12 Mi 



-2/3 



(5) 



We assume that the SMC is at a distance of 58.9 kpc 
( Kallivaya lil et al.|2006b| and has a stellar mass of Msmc* = 
3 x 10 s Mq ( Stanimirovic et al. 112004 1. The proper motion 



measurements of Kallivayalil et al. (2006b) give 

VsMC,tan 
VsMC.rad 



301 ± 52 kms" 
23±7kms" ] 



(6) 



From these values, the specific angular momentum of the 
SMC is 

~. RsMC VsMC,tan „.,/ Mvir \ ^ ,_s 

The angular momenta of the LMC and SMC are therefore 
strikingly close to each other. 

Piatek et al. (2008) have performed an independent 



analysis of the HST proper motion data and find results 
that are in general agreement with, but at the lower end 
of, the range found by K06. In particular, Piatek et al. 
find (Kan, Kad) of (346 ± 8. 5 km s" 1 , 93.2 ± S^kms" 1 ) and 
(259 ± 17kms~\ 6.8 ± 2.4kms _1 ) for the Large and Small 
MCs, respectively. Adopting these values changes the nor- 
maliza tion of Eqn. ([5]) to .52 and of Eqn. (7| to 0.46. Re- 



cently, Vieira et al. 



(20101 determined proper motions for 



the MCs using photographic and CCD observations from 
the Yale/San Juan Southern Proper Motion program span- 
ning a baseline of 40 years. They also confirm the K06 results 
and find proper motions for the SMC that are more consis- 
tent with the |Kallivayalil et ah] ( |2006b[ ) analysis than that 
of Piatek et al. 



2.3.2 Magellanic Clouds in the Millennium-II Simulation 

In order to identify z = analogs of the LMC, we first search 
our full MW sample for all subhalos that survive to z = 
and reside within a fixed fraction of R v i r of their host (see 
below). For each of these subhalos, we also find the mass 
of its main progenitor at all earlier times. To define analogs 
of the Magellanic Clouds in an iV-body simulation such as 
the MS-II, we need a way to link dark matter (sub)halos to 
galaxies. Recent work has shown that many observed sta- 
tistical properties of galaxies can be reproduced under the 
simple "abundance matching" assumption that stellar mass 
is a monotonically increasing function of the maximum dark 



matter mass a subhalo attains over its history (e.g., Conroy 



et al. I2006| | Wang et al.|2006| poster et al.|20To||Guo et al 



2010||Klypin et al.|201o| . Accordingly, we select as our fidu- 
cial sample the subhalo with the largest maximum main 
progenitor mass in each halo; we will refer to this sample as 
the first-ranked subhalos. The mass of this progenitor and 
the redshift at which it was reached are denoted M acc and 
Zucc', we will sometimes refer to these quantities as the in- 
fall mass and infall redshift. Note that subhalos selected in 
this manner are not necessarily the most massive subhalos 
in their hosts at z = 0. 

We identify LMC analogs within this sample by adding 
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two additional conditions. First, the host halos must obey 
M v i r G [1 — 3] x 10 12 Mq , corresponding to approximately 
35% of our full host halo sample; note that this is in the 
upper end of our host sample in terms of mass. Next, we 
use abundance matching as implemented by |Guo et al.| 
( |2010[ ) to estimate the maximum dark matter halo mass 
of the LMC. With our definition of halo mass, the in- 
fall halo mass corresponding to the LMC's stellar mass is 
M L mc(^cc) ~ 1.6 x lO n M . We define LMC analogs as 
first-ranked subhalos having infall masses within a factor of 
two of this mass, 8 < Af L Mc(^acc) [10 10 Mq] < 32, which 
allows for both scatter in the M* — Mhaio relation and un- 
certainty in the stellar mass of the LMC. Such subhalos 
are more massive than the typical most massive subhalo, 
in terms of maximum progenitor mass, of MW-mass halos 
(BK10; see Sj3]for a full analysis). 

To identify SMC analogs, we first find the surviving 
subhalos at z — having the second highest infall mass 
among all subhalos of their hosts; we will refer to this 
fiducial sample as the second-ranked subhalos. Using the 
SMC's stellar mass, we select second-ranked subhalos with 
masses between 4 x 10 10 and 1.6 x 10 11 Mq in hosts with 
M vir G [1 - 3] x 10 12 Mq as our SMC analogs. All SMC 
analogs therefore have at least 4000 particles at their maxi- 
mum, well above the ~ 1500 particle requirement that |Guo] 
et al.|(|2010l) show is necessary to adequately resolve the M acc 
halo + subhalo mass function at z = (see also |Wetzel fc] 
White|2010 |. 



0.0 0.2 0.4 0.7 1.1 1.9 4.0 



It is important to choose a suitable limiting radius for 
defining each halo's most (second-most) massive surviving 
subhalo. A natural choice is the redshift zero virial radius. 
We therefore consider the most (second-most) massive non- 
dominant subhalo within _R v i r at z — as our fiducial sam- 
ple. The precise choice of limiting radius used is not particu- 
larly crucial, so long as the radius is sufficiently large to cover 
much of the halo: we find similar results for 0.65 — 1.0 i? v i r - 

To summarize, we define two samples for both the LMC 
and the SMC: 
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• first-ranked subhalos (1st): 

the most massive subhalo, in terms of the maximum mass 
ever attained, located within _R v i r of a host's center at z = 0. 
This sample contains 7641 subhalos. 

• LMC analogs: 

the subset of 1st with 8 x 10 10 < M acc /Af Q < 3.2 x 10 11 
located in hosts with virial mass M vir G [1 - 3] x 1O 12 M . 
This sample contains 938 subhalos. 

• second- ranked subhalos (2nd): 

the second most massive subhalo, in terms of the maximum 
mass ever attained, located within R V [ T of a host's center at 
z = 0. This sample contains 7639 subhalos. 

• SMC analogs: 

the subset of 2nd with 4 x 10 10 < Af acc /M Q < 1.6 x 10 11 
located in hosts with virial mass M vir G [1 - 3] x 1O 12 M . 
This sample contains 840 subhalos. 



Figure 1. Top: Lookback time when the 1st ranked subhalo (solid 
histogram) or LMC analog (dashed histogram) first crossed the 
physical 2 = virial radius of its host. Vertical lines mark the 
medians of the distributions. Bottom: Cumulative version of the 
upper plot. 50% of LMC analogs were accreted by MW-mass hosts 
within the past 4 Gyr and 25% within the past 2 Gyr. 



3 ORBITAL PROPERTIES OF LMC 
CANDIDATES 

3.1 First crossing time 

The lookback time at which the LMC first crossed the phys- 
ical 2 = virial radius of the MW, moving inward - referred 
to hereafter as the first crossing time (ff c ) - serves as a useful 
discriminant between the early and late accretion scenarios 
laid out in B07. We therefore compute tf c for all first-ranked 
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Figure 2. Distribution of (specific) angular momentum for LMC 
analogs accreted early (more than 4 Gyr ago; magenta) and late 
(within the last 4 Gyr; black), as well as very early (more than 
8 Gyr ago; gray) and very late (less than 2 Gyr ago; cyan). Also 
plotted is the cumulative distribution for all lst-ranked subha- 
los, independent of accretion epoch (black dotted line). The gray 
shaded region corresponds to the observed angular momentum of 
the LMC, including ±2cr errors on the tangential velocity, for a 
Milky Way mass in the range [1 — 3] X 10 12 Mq . Late-accreted 
LMCs typically have higher specific angular momentum, due to 
both the inefficiency of dynamical friction over ~ 3 Gyr time- 
scales and the higher specific angular momentum of late-accreted 
subhalos. 



halos in our full MW sample to investigate whether there is 
a preferred accretion epoch for LMC-like objects. 

The distributions of tf c for all first-ranked subhalos and 
for those in the LMC analog subsample are shown in Fig. [T] 
as solid and dashed histograms, respectively. The full dis- 
tribution is bimodal, with a primary peak at t « 8 Gyr 
and a secondary peak at t = Gyr; the median value is 
tic = 5.8 Gyr. The distribution for the LMC-mass sample is 
markedly different. There is no prominent peak at t{ c ~ 8 
Gyr; instead, the distribution rises continuously from large 
lookback times (high redshift) to the present day, and the 
median of the distribution lies at tf c — 3.9 Gyr. The dif- 
ferences between the two distributions reflect the relatively 
high masses of the LMC analog sample: since these subha- 
los are more massive than the average first-ranked subhalo, 
they are accreted onto their host halos later, in the typical 
hierarchical manner expected in the ACDM cosmology. 

Candidates in the LMCs expected mass range are usu- 
ally accreted at fairly late times: only 12% have tf c > 7.5 
Gyr (2f c > 1), while approximately 30% have been accreted 
within the past 2 Gyr. Such numbers favor the late accretion 
scenario for the LMC. Even taking all first-ranked subhalos 
(the solid line in the lower panel of Fig. fIJ, we find that ~ 
70% were accreted since 2 = 1. First-ranked subhalos that 
could have completed several orbits are rare. 



1.0 




Figure 3. Distribution of orbital energies for the LMC analogs 
relative to the energy of a circular orbit at the host's virial radius. 
Blue and green shaded regions correspond to hosts of [1 — 2] X 10 12 
and [2 — 3] X 10 12 Mq, respectively. Unbound orbits (E < 0) are 
extremely rare, as are orbits consistent with a low-mass ( < 1.5 X 
10 12 M @ ) MW. 

3.2 Angular momentum 

Figure [2] shows the cumulative distribution of specific an- 
gular moments^] j = iiVtan, normalized by the virial value 
jvir = Rvii V v i v , for the full lst-ranked subhalo sample (dot- 
ted line). Also plotted is the same quantity for LMC analogs, 
split based on their accretion times: less than 2 and 4 Gyr 
ago (cyan and black lines, respectively), and more than 4 
and 8 Gyr ago (magenta and gray lines). The gray shaded 
region marks the range of j = j/j v ir allowed for the LMC 
based on a MW mass of 10 12 < M vil < 3 x 10 12 Mq and in- 
cluding ±2cr on the LMCs tangential velocity. This allowed 
range is fairly broad and shows that the angular momen- 
tum of the observed LMC is, in fact, fairly typical of LMC 
analogs in the MS-II. 

Comparing the angular momentum of LMC analogs 
based on accretion epoch, we find that LMCs accreted > 
4 Gyr ago have much lower angular momentum than those 
accreted at later times. In fact, 50% of those LMCs accreted 
>4 Gyr ago have angular momentum lower than the ob- 
served range for the LMC. On the other hand, half of the 
LMC analogs accreted within the past 2 Gyr have angu- 
lar momentum within the shaded region, supporting a late 
accretion scenario for the LMC. 

3 The angular momentum is computed with respect to the host's 
center, defined by the location of the gravitational potential min- 
imum. The host's velocity is determined by averaging over all 
particles in the host's main subhalo. We have also tried comput- 
ing the velocity with respect to only the most bound particles in 
the host subhalo or only those particles within 10 or 25 kpc and 
found that the results presented here are unchanged. 
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Early accreted (> 4 Gyr ago) LMC analogs in low-mass 
MWs are strongly disfavored: even in the most extreme sce- 
nario of M vir — 3 x 10 12 and Vtan,LMC ~ 331 kms^ 1 (2a 
lower than the central value of K06), we still find that over 
50% of LMC analogs have a lower value of j than for the 
actual LMC. A more realisitc scenario - M V - 1T — 1.6 x 10 12 
and Kan, lmc = 367 km s -1 - results 90% of early-accreted 
LMCs having lower angular momentum than is observed at 
z = 0. 

Regardless of accretion epoch, approximately 30-35% of 
LMC analogs fall in the gray shaded region: the LMCs an- 
gular momentum is not atypical in a cosmological context. 
This result dismisses previous assertions that tidal torques 
from M31 are needed to explain the orbital angular momen- 
tum of the LMC (e.g., |Raychaudhury fc Lynden-Bell||1989 



Shuter 1992 Byrd et al. 19941. Originally, concerns about 



the LMCs angular momentum arose because the orbital 
plane of the MCs is polar to the disk plane of the MW, while 
the LMC has orbital angular momentum that is at least as 
much as that of the MWs thin disk ( Fich fc Tremaine|19"91 



Lin et aL~| l995 Sawa & Fujimoto 2005). This is potentially 
difficult to explain in an early accretion scenario: torques 
from the MWs disk certainly could not explain the angular 
momentum of the LMC in this configuration, leading to the 
search for alternate potential perturbers. While our analysis 
does not extend to the likelihood of a polar orbital orienta- 
tion, we do find that the angular momentum of the LMCs 
orbit is more typical in a recent accretion scenario. Over 
such short interaction timescales, torques from the host are 
largely irrelevant and our result should hold regardless of 
the orientation of the orbit. 



3.3 Energy 

Figure [3] shows the orbital energy distribution of the MS- 
II LMC analogs, normalized by the energy of a circu- 
lar orbit at the host's virial radius. A striking feature of 
Fig. [3] is that nearly all subhalos are on bound orbits 
(E = E/ Scire (Rvii) > 0). The implications of the energies 
of subhalos' orbits for the LMC are highly sensitive to the 
virial mass of the MW: less than 10% of LMC analogs have 
orbits as energetic as that of the observed LMC if the MWs 
virial mass is smaller than 2 x 10 12 Mq (blue shaded re- 
gion), while a substantial fraction have orbits that match 
the energy of the LMCs orbit if the MWs halo lies between 
[2 - 3] x 10 12 M e ( green shaded region). Even if we take 
both radial and tangential velocities that are la (2a) lower 
than K06's mean values, we still find that only 2.3% (7.9%) 
of orbits are more energetic than observed for the LMC in 
a halo of 1.5 x 10 12 M Q . 

In Figure [4] we show the cumulative distribution of E 
for the LMC analog sample, split by accretion epoch (solid 
curves); the distribution for the full LMC analog sample 
(Fig. |3| is shown as the dotted curve. There is a marked 
difference when looking at early (>4 Gyr) versus late (< 4 
Gyr) accreted LMCs. Early-accreted LMCs tend to be on 
much more bound orbits (higher values of E), while late- 
accreted LMCs are less bound to their host halos (though 
virtually all are still bound, formally). 

Using the mean velocities from K06, we find that a 
Milky Way of mass (1, 2, 3) x 10 12 Mq corresponds to E = 




1 2 



Figure 4. Orbital energies for the LMC analog sample (dotted 
curve), as well as results split by tf c (solid curves). The early- 
accreted LMCs (> 4 Gyr, shown in magenta) are typically more 
bound to their host than are the late-accreted LMCs (< 4 Gyr, 
shown in black). This trend is even stronger for LMCs accreted 
over 8 Gyr ago (gray) versus those accreted within the past 2 Gyr 
(cyan). 



(-0.53, 1.21, 1.96). If the Milky Way's halo mass does not 
exceed 2 x 10 12 Mq, the early accretion scenario is strongly 
disfavored: there are vanishingly few early-accreted LMCs 
in the MS-II having E < 1.2. Fig. [I] also strongly disfa- 
vors any LMC accretion scenario for a 10 12 Mq MW (which 
puts the LMC on an unbound orbit). Approximately 25% of 
the late-accreted (and 7% of the early-accreted) LMCs have 
binding energies lower than the observed LMC for a MW 
mass of 2 x 10 12 M Q . The energetics favor a combination of 
a massive MW halo and a late-accreted LMC. 



3.4 Eccentricity 

To further explore the typical orbital properties of surviving 
LMC analogs as a function of accretion epoch, we consider 
the distribution of orbital eccentricities for these systems. 
Orbital eccentricity e is here defined as a combination of 
the pericenters r p and apocenters r a of orbits: 



r a + r v 



(8) 



(We assign unbound orbits an eccentricity larger than 1.) 
With this definition, e < 0.5 corresponds to fairly circular 
orbits (r a /r p < 3), while e = indicates a perfectly circular 
orbit. 

Figure[5]shows that early-accreted LMC analogs tend to 
be on orbits that are substantially more circular than those 
of late-accreted LMC analogs. Only 20% of LMCs accreted 
within the last 2 or 4 Gyr (cyan and black solid lines, respec- 
tively) have e < 0.5, a value that is met by approximately 
50% of LMCs accreted over 4 Gyr ago and by 60% accreted 
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Figure 5. Orbital eccentricity for all LMC analogs (dotted line), 
as well as split by accretion epoch (solid lines). Early-accreted 
LMCs (gray and magenta lines) have orbits that are significantly 
more circular than do late-accreted LMCs (cyan and black). A 
MW of mass (1.5,2,3) X 10 12 M corresponds to an LMC ec- 
centricity of (0.92, 0.79, 0.62), while a 10 12 Mq MW results in an 
unbound orbit (e > 1). 



more than 8 Gyr in the past. Although late accretion does 
not a priori mean that the Clouds cannot have completed 
multiple pericentric passages, these data on eccentricities 
provide further constraints. The mean eccentricity for the 
late-accreted LMC analogs is ~ 0.7. Given that the percen- 
ter of the LMCs orbit is ~ 45 kpc, the resulting apocenter is 
~ 260 kpc. The MCs clearly could not have completed mul- 
tiple pericentric passages on such an orbit within the past 4 
Gyr. We therefore conclude that the MCs were most likely 
accreted within the past 4 Gyr and are on their first passage 
about the MW 

Fig. [5] reinforces how unlikely it is to find an object 
with the LMCs infall mass and orbit in a 10 12 Mq halo at 
z = 0, regardless of accretion epoch, as the LMC has e > 1 
according to our definition - i.e., the LMCs orbit is unbound 
- for this mass. If the MW does have such a low-mass halo, 
then the LMC was certainly accreted within the past 4 Gyr. 
Moreover, less than 10% of orbits - independent of accretion 
epoch - have e > 0.92, which corresponds to a halo mass of 
< 1.5 x 10 12 Mq. It is therefore quite unlikely that the Milky 
Way has a mass of less than 1.5 x 10 12 Mq. (For reference, a 
Milky Way halo of mass 2x 10 12 Mq results in an eccentricity 
of 0.79 for the LMC, while a mass of 3 x 10 12 Mq gives an 
eccentricity of 0.62.) 



4 FREQUENCY OF LMC/SMC ANALOGS 
ABOUT MW HOSTS 

In this section, we examine the mass distribution of the sam- 
ples defined in § |2.3.2| Our goal is to determine how typical 



the infall masses of the MCs are relative to the full first- 
and second-ranked subhalo samples and as a function of the 
host mass. 

The upper left panel of Fig. [6] shows the distribution 
of masses for first-ranked subhalos (solid histogram) and 
LMC analogs (dashed histogram), relative to the redshift 
zero virial mass of their hosts. The distribution of masses 
for first-ranked subhalos is fairly broad and peaks at Mist ~ 
0.03Af v i r (z = 0). The mass distribution of LMC analogs is 
substantially narrower and peaks at a much higher mass, 
Mlmc ~ 0.1M V ii(3 = 0). The difference between the two 
distributions indicates that the LMC is more massive than 
the typical first-ranked satellite galaxy in a Milky Way-mass 
halo, a conclusion also reached in BK10. 

A similar situation exists for the SMC, which is shown 
in the upper right panel of Fig. [6] While the distribution 
of /U2nd peaks at ~ 0.01, the distribution of /xsmc peaks at 
~ 0.03. This shows that the SMC is also more massive than 
the second-ranked galaxy in a typical MW-mass halo. 

The solid line in the lower left panel of Fig. [6] shows the 
cumulative version of Mi st (z a cc) / M v [ T (z = 0). The combined 
shaded region encompasses the full LMC analog sample; the 
range of M a cc/Mvir(-Z — 0) corresponding to M V \ T = 10 12 Mq 
is shown in blue, while the range for 3 x 10 12 Mq in shown in 
gray. The lower right panel of Fig.[6]plots the same quantities 
for our sample of second-ranked subhalos and SMC analogs. 

Recall that there are 2658 host halos in the mass range 
M vir 6 [1 - 3] x 10 12 Mq . If we define MC analogs strictly 
in terms of mass, we thus conclude that approximately 35% 
of MW-mass halos host an LMC analog and 32% host an 
SMC analog within R v i r . These numbers are for a specific 
range of MW halo masses, however, and are sensitive to 
the precise mass of the MW (BK10). A halo of 10 12 Mq 
has a w 20% chance of hosting an LMC analog, and less 
than a 10% chance to host an SMC analog. A 3 x 10 12 Mq 
MW makes L/SMC analogs much more common, as approx- 
imately 40% of such hosts have L/SMC analogs. [Note that 
these numbers are likely to be upper estimates, as we have 
used very conservative estimates on the errors for the map- 
ping between M + and M acc .] 

In a search for LMC analogs in the seventh data release 
of the Sloan Digital Sky Survey ( York et al.|2 000 Abazajian 



et al.| [2009), Tollerud et al. (in prep.) determine that ap- 



proximately 40% of isolated hosts with luminosities similar 
to that of the MW also have an LMC analog (r-band magni- 
tude between -17.5 and -20) located within a projected sepa- 
ration of 250 kpc. Although the selection criteria of Tollerud 
et al. differ from those used here, this result appears consis- 
tent with our findings. 

The mass of the first-ranked subhalo correlates strongly 
with its first crossing redshift Zf c , which is shown in Fig. [7] 
The median (middle curve) decreases from zt c ~ 1.4 at /1 = 
0.01 to 0.4 at /1 = 0.1 (70% of the distribution is contained 
within the dashed curves). The shaded region in Fig.[7]shows 
the allowed range for the LMC when using the abundance 
matching assumption (see Fig. |6| 

An alternate way of looking at the dependence of /x(lst) 
on Zf c is to divide the 1st ranked subhalos into two samples: 
those accreted early, defined here as zf c > 0.4 (or equiva- 
lently, a look-back time of t > 4 Gyr, corresponding to the 
local minimum in Fig. [I]), and those accreted at late times 
(Zf c < 0.4 or t < 4 Gyr). The result of this split is shown in 
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Figure 6. The distribution of masses of first-ranked subhalos (left panels) and second-ranked subhalos (right panels) with respect to 
the virial mass of their hosts. Upper panels: Distribution of Mi st /M v i r (z = 0) and AfLMc/-^vir(^ = 0) (solid and dashed histograms, 
left panel), and of M2 n d/M v i r (z = 0) and MgMC / Mvir{z = 0) (solid and dashed histograms, right panel). Lower panels: Cumulative 
distributions of solid curves from upper panels. The shaded regions correspond to the observed range of infall masses for the LMC (left 
panel) and SMC (right panel) assuming a MW mass of either 10 12 Mq (blue region) or 3 X 10 12 Mq (gray region). 



Fig. [8] and reinforces the result of Fig. [7| 1st ranked subha- 
los that have joined their z = host halos within the last 4 
Gyr tend to be a factor of ~ 3 more massive at infall than 
those that joined their host halo more than 4 Gyr ago. Since 
the LMC has been shown to be more massive than the typ- 
ical 1st ranked halo (see § H, it is more likely to have been 
accreted at late times. 

In order to assess whether a late or early accretion sce- 
nario is more plausible for the LMC, we need to also account 
for the survivability/mass loss expected for LMC analogs, 
which will depend on the accretion epoch and the specific 
orbit of the subhalo. Figure [9] shows the distribution of ac- 
cretion mass relative to redshift zero dark matter masflfor 
LMC analogs. Results are plotted for early (more than 4 
Gyr ago; cyan) and late (less than 4 Gyr ago; black) accre- 
tion epochs, as well as for very early (more than 8 Gyr ago; 
gray) and very recent (less than 2 Gyr ago; cyan) accretion 



This is the bound mass determined by the SUBFIND algorithm. 



epochs. There is a pronounced, and not unexpected, trend 
for stronger mass loss in earlier-accreted LMC analogs. The 
most recently accreted LMCs are typically a factor of 1.5-4 
less massive at z = than at accretion (60% confidence in- 
terval), whereas the earliest accreted LMCs tend to be 3.5 
to 30 times less massive. It is unlikely that the LMC could 
have undergone strong tidal stripping without losing much 
of its gas, which is at odds with observations. 

Massive satellites typically do not survive for long 
before merging with their hosts: the dynamical friction 
timescale for a 1:10 object at z = 1 is approximately 5 
Gyr ( Boylan-Kolchin et al.|2008 1 , shorter than the time be- 
tween z = 1 and the present day. This is corroborated by 
Stewart et al. (20081, who showed that LMC-mass objects 
(M ~ 10 11 Mq/i -1 ) typically do not survive for more than 
~3 Gyr after accretion (their figure 5). BK10 examined the 
accretion epochs of massive subhalos in MW-mass hosts and 
found the same trend (their Fig. 12). Both groups argued, as 
we do here, that this lends support to a first passage scenario 
for the LMC. 
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Figure 7. Redshift Z{ c where the first-ranked subhalo first crossed 
the z = virial radius of its host as a function of M\ at /M v - lv (z = 
0). The middle line shows the median, while 70% of the distribu- 
tion is contained between the upper and lower lines. The shaded 
region(s) are the same as in Fig. [6] 



Figure 9. The cumulative distributions of infall mass relative 
to z = mass (in dark matter) for LMC analogs. Each curve 
corresponds to LMC analogs with different accretion epochs. The 
most recently accreted candidates (cyan line) have typically lost 
of order half of their M acc , whereas the earliest accreted LMCs 
(gray line) have lost almost 90% of their dark matter mass. 
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Figure 8. Distribution of masses for 1st ranked subhalos accreted 
early (more than 4 Gyr ago; black) and late (within the last 4 Gyr; 
magenta). Note that each distribution is separately normalized to 
unity; the number of first-ranked satellites with early accretion 
times is approximately 2.4 times greater than those with late 
accretion times. Those accreted at late times tend to be a factor 
of ~3 more massive at infall than those accreted at early times. 



5 LMC-SMC PAIRS 

The small separation in position (~20 kpc) and velocity 
space (~100 km/s) between the two Magellanic Clouds 
makes it unlikely for them to be coincidental neighbors. Fur- 
thermore, the existence of (1) a bridge of gas connecting the 
two galaxies and (2) the Magellanic Stream, a stream of HI 
gas that trails behind the MCs over 150° across the sky, also 
strongly suggests that they have interacted for at least some 
time in the past (Besla et al. 2010). 

These observations lead to a number of questions re- 
garding the MC system, including: What is the probability 
that the MCs were accreted together and survive as a bi- 
nary today? Is there a preferred accretion epoch or mass 
ratio for such a pair? How likely is it that the Clouds are 
only an apparent binary system today rather than a true 
binary system? 

To this end, we consider the difference in accretion 
epoch between each LMC analog and the corresponding 
second-ranked halo about the same host. Figure [10] shows 
the median difference (solid line) in tf c for the LMC analogs 
and second-ranked subhalos as a function of U c for the LMC 
analog, as well as the 10% (dotted) and 25% (dashed) quan- 
tiles. It is indeed possible to find first and second-ranked sub- 
halos that are accreted as pairs (within 1 Gyr of each other). 
This probability depends on the LMC accretion epoch, how- 
ever, and does not provide information about whether the 
accreted pairs can remain as a binary to the present day. 
The number of pairs accreted > 4 Gyr is approximately 
three times larger than the number accreted within the last 
4 Gyr. This is likely a result of 1 Gyr being a much larger 
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Figure 10. Difference in first crossing times between LMC 
analogs and the second ranked subhalo in each system as a func- 
tion of the first crossing time of the LMC. (All times are lookback 
times.) The dotted (dashed) lines show the 10-90 (25-75) per- 
centiles of the distribution, while the solid line shows the median. 
For LMCs with t fc ~ 4 Gyr, almost 40% of SMCs were accreted 
within ±2 Gyr of the LMC, and 25% within ±1 Gyr. 



fraction of a typical orbital time at higher redshifts than 
today. 

We can also compute how likely it is to find MC-like 
objects around MW-mass hosts in a binary system at z = 0. 
To do so, we calculate the separation in position and velocity 
space between the LMC analog sample and the correspond- 
ing second-ranked subhalo for the same host. Systems having 
|Av| < 150km/s and \AR\ < 50kpc/h are considered plau- 
sible binaries. Table 1 lists a number of properties for the 
23 identified LMC/SMC analogs, including the first cross- 
ing time for each Cloud. Recall that the LMC analog sample 
contains 938 halos; it is thus possible, though not probable 
(~ 2.5%), that an LMC analog and a second ranked halo be 
found in a binary system about a MW-mass host today. 

Since SMC analogs represent a subset of the second 
ranked subhalo sample, these numbers also indicate that 
LMC / SMC analogs accreted at similar epochs are not likely 
to exist as binaries at the present day. The binaries that do 
survive to z = tend to have the SMC in an eccentric orbit 
about the LMC. This is in agreement with the proper motion 



measurements of the SMC by Kallivayalil et al. ( 2006b I and 



Vieira et al. ( 2010 1, who find a high relative velocity between 
the Clouds. An eccentric SMC orbit about the LMC is also 
required in the Magellanic Stream model proposed by |Besla| 
et al. ( 20101). The present study shows that such an orbital 



configuration is cosmologically expected. 

A few specific cases from Table [l] stand out in particu- 
lar. The first three rows have the pairs that were accreted 
most recently. All three of these systems have more angular 
momentum than the true LMC, highlighting that there may 
be no issue with the large magnitude of the LMCs angu- 



Table 1. Properties of LMC-SMC binaries. Column 1: ratio 
of masses of most massive progenitors; column 2: mass ratio at 
z = 0; column 3: first crossing time for LMC; column 4: first 
crossing time for SMC; column 5: separation between LMC and 
center of host; column 6: 3-D velocity of the LMC, relative to 
host. 
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a Note that the calculation of subhalo masses at redshift zero is 
sensitive to the location of the subhalos within their hosts. 



lar momentum if it was accreted recently. Additionally, all 
are on fairly energetic orbits: with E/E V - 1T w 1, these are 
the most energetic among all 23 binary candidates. In fact, 
most of the other candidates have E/E v { r > 2, which places 
them on orbits that are improbably bound (relative to the 
observed LMC) even for a MW of 10 12 , M© (see Fig. [|]). 
With the exception of the oldest binary candidate (the final 
row of Table IT]), all systems with tf c (LMC) > 5 Gyr have 
low angular momentum, likely due to losses via dynamical 
friction. Although the volume of the MS-II does not provide 
a vast sample of possible binaries, those that do exist with 
orbital energetics similar to that observed for the Clouds are 
accreted at late times. Specifically, the systems highlighted 
in the first three rows of Table 1 illustrate that it is possible 
for a binary LMC/SMC to be accreted very recently on a 
high angular momentum/energy orbit 

A further point of interest is that the pairs with first 
crossing redshifts that differ by > 2.5 Gyr between the LMC 
and the corresponding second ranked subhalo are all cases 
where the LMC is within 100 kpc of the host. If these sys- 
tems were true binaries, then there should also exist exam- 
ples where the two subhalos had discrepant accretion epochs 
and are located at large distances from the host today. Since 
no such examples exist, it is likely that these are chance asso- 
ciations of two satellites near pericenter than true binaries. 



© 2010 RAS, MNRAS 000,[TJ{T4] 



Magellanic Clouds in a ACDM Universe 11 



1.0 I 1 1 1 1 1 I 1 1 1 1 r 

0.20 - 




0.1 1 1 1 1 0.00 I 1 1 1 1 ' " ' I 1 

0.5 1.0 1.5 2.0 2 4 6 8 10 12 

1 +<5 5 lookback time [Gyr] from z=0.21 



Figure 11. Dependence of Zf c on the large-scale overdensity, mea- 
sured via gaussian smoothing on 5 /i -1 Mpc spheres, for the LMC 
analogs. There is no obvious dependence of Zf c on the large-scale 
overdensity, indicating that the large-scale environments of MW- 
mass hosts does not strongly influence the accretion epochs of 
their most massive subhalos. 



Figure 12. Distribution of first crossing times for the sample of 
MW-mass hosts selected at z = 0.21 (analogous to upper panel of 
Fig. Q . This distribution is similar to the distribution for z = 
hosts, suggesting that changing as from 0.9 to 0.815 should not 
affect our conclusions on the probability of the LMC being on its 
first pass about the Milky Way. 



6 DISCUSSION 

Using the Millennium-II Simulation, we have investigated 
ACDM predictions for orbital properties and accretion his- 
tories of the MCs in the context of the updated proper mo- 
tion measurements of Kallivayalil et al. (2006a, b). In this 
section, we further explore how environment or cosmologi- 
cal parameters may influence our results, the likelihood that 
the MCs are on their first passage about the MW, and the 
masses of halos hosting subhalos with LMC-like separations 
and velocities. 



6.1 Milky Way sample 

It is important to investigate whether the trends shown in 
the previous sections have any systematic dependence on 
any properties of the host halos. In particular, halos resid- 
ing in low density environments have different accretion his- 



tories than those in high density regions (Gottlober et al 



2001 Maulbetsch et al. 2007 Fakhouri & Ma 2010 1, an effect 



that could potentially bias our results on accretion epochs 
of LMC analogs. We therefore plot how the large-scale en- 
vironment of a halo influences the first crossing redshift of 
LMC analogs in Fig. |11| This plot shows that Zf c is essen- 
tially independent of environment as measured by the dark 
matter overdensity, smoothed with a Gaussian filter of width 
5 Mpc. Our results should therefore be insensitive to the 
environment of the host halo. 

The typical accretion epoch of LMCs could also be af- 
fected by the choice of cosmology in the MS-II, which has 
as that is approximately 10% higher than the current best- 
fitting value of « 0.81 (Komatsu et al. 2009[ ). [This difference 



is minor in terms of the number of MW-mass halos found 
at z = 0, as it affects the abundances of such halos by 10%.] 
While a quantitative understanding of the effects of varying 
ag would require running a completely new simulation, we 
can easily estimate the qualitative effect by noting that low- 
ering as results in later formation of halos of a given mass. 
We therefore expect that any changes due to reducing as 
would tend towards even later accretion epochs for LMC 
(and SMC) analogs. 

Alternatively, we can note that cr g (z = 0.21) « 0.81 
for the MS-II cosmology; the difference between the distri- 
bution of if c for halos defined at this epoch and at z = 
should therefore inform us about potential differences due 
to changes in the power spectrum normalization. The dis- 
tribution of first crossing times for lst-ranked subhalos and 
LMC analogs of halos selected from the MS-II at z = 0.21 
(using the same criteria described in §|2| is shown in Fig. 12 
Comparing this distribution with the distribution of tf c in 
halos selected at z = (the upper panel of Fig.[T|, we can see 
that little changes between z — and z = 0.21. Accordingly, 
reducing as from the Millennium and MS-II value of 0.9 is 
not likely to strongly affect our findings on the accretion 
epochs of LMC-like satellites. 



6.2 Are the Magellanic Clouds on their first 
passage about the Milky Way? 

A number of lines of observational evidence support the idea 
that the Magellanic Clouds are making their first pericentric 
pass about the Milky Way. This scenario explains why two 
gas-rich satellites reside at small galactocentric distances - 
similar satellites are typically found at much larger distances 
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from the MW or M31 (van den Bergh|2006[ ) - as the Clouds 
would not have had sufficient interaction with the MW to 
have lost their gas by some combination of tidal stripping, 
harassment, and ram pressure stripping ( Mayer et al.|2 006) 



1.0 



Similarly, the unusually blue color of the LMC (James & 
|rvory|[2~010| Tollerud et al., in preparation) means that it 
must have retained a substantial amount of star-forming gas, 
which is difficult to understand if the MCs have completed 
multiple orbits about the MW. The existence of stellar pop- 
ulations extending as far as « 20 kpc from the LMC's dy- 



namical center (Munoz et al. 2006 Majewski et al. 2009 



Saha et al.||2010[ ) is also an indication that the LMC has 
not interacted strongly with the MW. Finally, |Besla et al.| 
(20101 have shown that the Magellanic Stream may origi- 



nate from a tidal interaction between the MCs themselves, 
a model that requires the MCs to be a recently-accreted 
binary system. 

B07 first examined the possibility that the MCs have 
been recently accreted by the MW using an orbital analysis 
constrained by the new HST proper motion measurements. 
Uncertainties in modeling the MW meant that a scenario 
in which the Clouds were accreted at early times cannot 



be ruled out by such an analysis (e.g., Gardiner & Noguchi 



p96l[BesIa et al.||2007| |Shattow fc Loeb||2009[ ), and refine 
ments in the error space of the proper motions are unlikely 
to improve this situation. (We note, however, that early ac- 
cretion models assume a static MW halo potential over a 
Hubble time, which is unrealistic.) We have computed the 
likelihood of a first passage scenario using a large sample 
of MW-mass halos from a high resolution cosmological sim- 
ulation of the ACDM cosmology. Our results put the first 
passage scenario on even firmer ground. We have showed 
that it is highly improbable for surviving satellites with in- 
fall masses similar to that of the LMC to have been accreted 
at z > 1, rendering an early infall scenario unlikely based 
on mass considerations alone. We further found that 25% of 
surviving LMC analogs have been accreted within the past 
2 Gyr and that the energetics of the LMC orbit are strongly 
inconsistent with the properties of LMC analogs accreted > 
4 Gyr ago. Finally, we have showed that recently-accreted 
LMCs are incapable of making multiple pericentric passages 
by the present. Taken together, these results demonstrate 
that it is quite likely that the MCs have recently joined the 
MW and are currently making their first pericentric pass; a 
very recent accretion ( < 2 Gyr ago) is also cosmologically 
plausible. 

6.3 Mass of the Milky Way 

Uncertainties in the mass of the Milky Way's halo play a sub- 
stantial role in placing the LMC and its orbit in a cosmolog- 
ical context, as has been shown repeatedly in our analysis. 
A low mass halo (M vir « 1O 12 M ) means that the LMC 
is fairly unusual in terms of its (high) mass and that it is 
very unusual in terms of its (energetic) orbit. Both the mass 
and orbital energy of the LMC are more typical for halos of 
M vir >2 x 1O 12 M . 

We can also take an "inverse" view and ask, in what 
mass dark matter halo do objects with masses, velocities, 
and halo-centric distances similar to the LMC reside? To this 
end, we build a sample of LMC analogs with no constraint 
on the properties of the host halo (note that this differs 



M acc G [8-32] xlO 10 M s 
35<i?<65 kpc 
300<V<420 km/s 

2<0.3 (495 subhalos) 



V 
> 




Figure 13. Distribution of halo masses in the MS-II having a 
satellite similar to the LMC: a mass of 8 x 10 10 < Af acc /Af < 
3.2 X 10 11 M , at a separation of 35 < R < 65 kpc from its 
host, and with a total velocity of 300 < V < 420kms~ 1 . Of 
all the hosts with such a subhalo at z < 0.3, over 90% have 
M vir > 2 x 10 12 M Q . 



from the mass-selected samples used up to this point). We 
merely require the following properties of the subhalo: (1) 
0.8 < M acc [10 11 M@] < 3.2, (2) 35 < R < 65 kpc from its 
host, and (3) 300 < Vtot < 420kms _1 . Since these criteria 
are somewhat restrictive, we search for hosts in all 10 MS- 
II snapshots with z < 0.3. We find 495 subhalos matching 
our search criteria; the distribution of host halo masses for 
these matches are shown in Fig. |13| The figure confirms that 
satellites with properties similar to those of the LMC are 
unlikely to reside in host halos with M vlr < 1.5 x 10 12 Mq. If 
we further require that M V - 1T < 3 x 10 12 M , we still find that 
approximately 70% of hosts of LMC-like subhalos reside in 
hosts with Mvir > 2 x 10 12 M©. 



7 CONCLUSIONS 



The new HST proper motions for the MCs ([Kallivayali l et ah] 
2006a|b I have forced us to re-evaluate our understanding of 
their orbital history about the MW. The canonical picture, 
wherein the MCs are on a quasi-periodic, slowly decaying or- 
bit around the MW, has been thrown out. We are left instead 
with two possibilities: (1) The MCs are on their first pas- 
sage about the MW (late accretion); or (2) the MCs joined 
the MW > 8 Gyr ago and are now on a highly eccentric 
orbit, having already completed at least one passage about 
the MW since infall (early accretion). 

In this work, we have addressed the likelihood of these 
two scenarios by using the MS-II to place the MCs in a cos- 
mological context in terms of their accretion epoch, orbital 
properties, and masses. Our primary results can be summa- 
rized as follows: 
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• LMC analogs are accreted preferentially at late times: 
only 15% have tf c > 7.5 Gyr (zf c > 1), while approximately 
30% have been accreted within the past 2 Gyr. Such numbers 
favor the late accretion scenario for the LMC. 

• The LMC's angular momentum is not anomalously 
high: 30-35% of all LMC analogs have specific angular mo- 
mentum matching that of the real LMC if the mass of the 
MW lies between [1 — 3] x 10 12 Mq. The angular momentum 
of the LMC is more typical of halos at the massive end of 
this range than of those at the low-mass end. 

• It is exceedingly unlikely for the LMC to be on an un- 
bound orbit: If the mass of the MW is less than 2 x 10 12 Mq , 
the LMC has an orbit that is more energetic than 90% of 
comparable systems (adopting the mean radial and tangen- 
tial velocities of K06). 40% of MW systems with M v ; r € 
[2 — 3] x 10 12 Mq have an LMC analog with orbital energy 
comparable to that of the LMC. It is highly unlikely for 
LMC-like subhalos to be on unbound orbits, which is the 
case for a low-mass MW, and none of the early-accreted 
LMCs are on unbound orbits. The conclusion that the LMC 
is in fact bound to a massive MW, and yet accreted re- 
cently, cautions against the use of backward orbital integra- 
tion schemes to determine the orbital histories of satellites 
over cosmic time. 

• Energetically, it is difficult to accommodate a scenario 
where the MCs have made multiple pericentric passages: 
LMCs accreted at early times are on mostly circular orbits, 
at odds with observations. LMCs accreted recently have not 
had time to complete more than one pericentric passage. 

• LMC and SMC-mass objects are not particularly un- 
common in MW-mass halos: In a refinement of the results 



presented in Boylan-Kolchin et al. (20101, we find that 20- 



32% of MW-mass halos host an LMC analog and 10-25% 
host an SMC analog. These results are consistent with the 
analysis of LMC analogs about MW type hosts located in 
the SDSS DR7 by Tollerud et al. (in prep.). The MCs be- 
come less typical if the host halo is lower in mass. 

• It is possible, but not probable, to find LMC-SMC bi- 
naries at z—0: We find a small number of Milky Way-mass 
systems (~ 2.5%) with LMC analogs have apparent LMC- 
SMC binaries. 

• Subhalos with properties similar to that of the LMC re- 
side preferentially in massive host halos: Out of all dark 
matter halos (without restriction on mass) hosting objects 
with masses, velocities, and separations similar to the LMC, 
only 10% have M vir < 2 x 10 12 Mq, and less than 5% have 
M vi[ < 1.5 x 1O 12 M . 

Overall, our results support a scenario in which the 
LMC is a recent addition (in the last 4 Gyr) to a fairly 
massive (M vir > 1.5 x 10 12 Mq) Milky Way. 
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